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ABSTRACT 

Phenological data as reported at the crop reporting district (CRD) level, 
and environmental data from the National Climatological Center appear to be 
viable sources of information for use in development and testing of an 
adjustable crop calendar model for winter wheat. These data were utilized 
for this study, in which generalized least squares techniques were applied 
for parameter estimation of functions to predict winter wheat phenological 
stage, with environmental values as independent variables. The independent 
variable investigated included daily maximum temperature (T^) , daily minimum 
temperature (T^) , daily daylength (Dj ) , and daily precipitation (P^). 

After parameter estimation, tests were conducted with variances from 
the fits, and on independent data. From these tests, it may generally be 
concluded that exponential functions have little advantage over polynomials. 
Precipitation was net found to significantly affect the fits. The Robertson’s 
"triquadratic" form, in general use for spring wheat, was found to show 
promise for winter wheat, but special techniques and care are required for 
its use. In most instances, equations with nonlinear effects were found to 
yield erratic results when utilized with daily environmental values as 
independent variables. Thus, as of this writing, the linear function of the 
form 



R * H. + H^T + H-T + H.D. 
1 2 x 3 m A L 


is recommended. Specific coefficients recommended for inclusion and testing 
in the LACIE project are: 


ii 


til 


Stage ... 
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Plant-Fme rge 

-0.014919 

0.0038970 

0 . 

0 . 

Erne rge-.To int 

-0.00039918 

0.00043509 

0 . 

0 . 

Joint -Head 

-0,216419 

0 . 

0 . 

0.019021 

Head-Soft Dough 

0.314583 

0 . 

0 . 

-0.018610 

Soft Dough-Ripe 

0.24471 l 

0.00462 l l 

0.0015439 

-0.022684 


Specific reeommendat tons tor further work, to be conducted during the 
second phase of the contract, include preparation and inclusion of further 
data in the least squares programs; preparation of more extensive testing 
programs and data, to Include investigation of the effects of using 
averaged environmental data for predictions; further work on the Robertson's 
”t r iquadrat ic" model; and variance propagation studies. 
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INTRODUCTION 

The Large Area Crop Inventory Experiment (LACIF) utilises winter wheat 
crop calendar, or biometeorological time scale, information in several ways. 
Among these, analyst interpreters use the information in training field 
selection and other researchers use the information in yield modeling efforts. 

To date the needed winter wheat crop calendar information has been pro- 
vided in one of two ways. In the first, historical records are used to com- 
pute an average date at which the crop will reach a specific stage. A 
"spread", computed statistically from the data, is included to account for 
yearly variations. In the second, functions found by Robertson (1)* are 
modified for use on winter wheat. These modifications are necessary be- 
cause the Robertson’s model was generated for spring wheat in the Canadian 
prairie provinces. The most commonly used procedure used to modify the 
Robertson’s model has been to calculate "multipliers" for winter wheat. 

Both of these techniques suffer from serious limitations. In the 
first, variations in planting dates, and the large variations inherent when 
averaging over many years, means that considerable uncertainties result in 
the estimated dates. For the second, it must be emphasized that Robertson’s 
model was generated for spring wheat, whose growth environment differs 
considerably from that of winter wheat. Spring wheat temperatures in gen- 
eral are on the rise throughout the growth and ripening of the crop. 

Day lengths in general increase during emergence and subsequently decrease. 
Winter wheat, on the other hand, is planted and emerges during a period of 


^Numbers in parentheses throughout this report refer to the reference list 
in the back. 
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declining temperatures and day lengths, undergoes a long winter dormancy 
period, and then completes its growth under conditions of increasing 
temperature and davlengths after spring grcenup. 

For these reasons it was felt desirable to seek an acceptable crop 
calendar for winter wheat specifically, based upon phenological reports, in 
which weather and environmental variables could be used to predict winter 
wheat growth stages. This was the overall goal of the research conducted. 
With this goal in mind the following research objectives were formulated: 

a. to define data sources which may be utilized in the formula- 
tion of such a model; 

b. to generate computer programs which may be used to investi- 
gate various mathematical models; 

c. to identify functional forms which are appropriate for des- 
cribing winter wheat growth stages; 

d. to obtain reliable estimates of parameters for these func- 
tions, and variance estimates of these parameters; 

e. to identify and investigate auxilliary problem areas, such 
as winter dormancy criteria; 

f. to provide recommendations for further studies. 


FUNCTION ESTIMATION 

The objective of the modeling process was to obtain a relationship 
predicting phenological stage numbers as a function of environmental 
variables. The phenological stages and the associated stage numbers used 
in this study are as follows: 
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Sta ge Number 

Planting 0 

Emergence 1 

Jointing 2 

Heading 3 

Soft-Dough 4 

Ripe 5 

The environmental variables selected were the same as those used by 
Robertson (1). These variables are daily maximum temperature, daily 
minimum temperature, and daylength, used to represent photoperiod. In addi- 
tion, precipitation was investigated as a possible independent variable, in an 
attempt to assess the possible importance of moisture in such a model. 

Data Sources 

Phenological data at the Crop Reuorting District (CRD) level were used 
in this study. Depicted in Figure 1 are the crop reporting districts used in 
the United States. Historical data from some Great Plains, Midwest, and 
Rocky Mountain states were obtained from NASA’s Johnson Space Center, Earth 
Observations Division. Table 1 shows a summary of the phenological data 
obtained which was available for the investigation. The dates used for 
modeling were the 50X, or median dates, converted to Julian Day Numbers, at 
which the crop reached each stage. Table 2 shows a typical report of the dates 
at which the crop reached each stage for a single year in one state. It 
should be noted that, throughout the United States, there exists a great 
variability as to winter wheat phenology reporting. As may be seen from 
Table 1, few states report all the phenological stages, and several report 
stages which differ from the standard stages. For example, Oklahoma reports 
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Figure I: Crop Reporting Districts in the United States 
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Table 1 

Summary of Pheno logical Data Available 


Plant -Emerge : 


Emerge-Joint : 


Joint-Head: 


Oklahoma^ - 19 6U , 65, 66, 6T, 68, 69 , TO, 71, 72, 73, 7^ 
Colorado - 1972, 73, 7 1 * , 75 
Idaho - 1973, 7 1 *, 75 

Oklahoma - 196U , 65 , 66, 67, 68, 69 , 70, 71, 72, 73, 7^ 

Colorado - 1972, 73, lb 

North Dakota 2 - 1970, 71, 72, 73, 7^, 75 

Montana 2 ’ 3 - 1971, 72, 73, 7**, 75 

Idaho 3 - 1973, 7*», 75 

Oklahoma - 1964, 65 , 66, 67 , 68, 69 , 70, 71, 72, 73, 7 1 * 
Colorado - 1971, 72, 73, 74 

Kansas - 1964, 65 , 66, 67 , 68, 69 , 70, 71, 72, 73, 74 

Idaho 3 - 1973, 74, 75 

Montana 3 - 1971, 72, 73, 74 

North Dakota - 1970, 71, 72, 73, 74, 75 


Head-Soft Dough: Oklahoma - 1964 , 65, 66, 67 , 68, 69 , 70, 71, 72, 73, 74 

Kansas - 1964, 65 , 66, 67 , 68, 69 , 70, 71, 72, 73, 74 

Idaho - 1973, 74, 75 

Montana - 1971, 72, 73, 74, 75 

Missouri - 1970, 71, 72, 73 


Soft Dough-Pipe: Oklahoma - 1964, 65 , 66, 67 , 63, 69 , 70, 71, 72, 73, 74 

Kansas - 1964, 65 , 66, 67 , 68, 69 , 70, 71, 72, 73, 74 

Idaho - 1973, 74, 75 

Montana - 1971, 72, 73, 74, 75 


Notes : 1.) "Acceptable Stand" reported. Emergence estimated by 

interpolation. 

2. ) Emergence not reported. Spring greenup period only used in 

regression. 

3. ) "Boot" reported. Jointing estimated by interpolation. 


Table 2: Typical Report of Winter Wheat Phenology for CRD's 


S’ a .i . _ KANSAS 
C o»# V2.»r: 72 
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"acceptable stand" rather than emergence, and both Idaho and Montana report 
"boot" rather than jointing. 

As independent variables, values of daily maximum and minimum temperature, 
daily precipitation, and daylength may be used. Values for temperatures and 
precipitation are available on a daily basis from the National Climatological 
Center, NOAA, Asheville, North Carolina. As listed in Table 1, records were 
obtained for the states and years for which phenological data were available. 
Tables 3 and 4 show typical listings of temperatures and precipitation, 
respectively, as reported by the National Climatological Center. The value 
of daylength may be calculated as a function of latitude and Julian Day 
Number. The empirical equation obtained by Stuff (2) was used for this pur- 
pose in this investigation in which a tangent function is utilized for 
latitudes less than 40°, and a tangent-squared function for latitudes above 
40°. 

Functional Forms 

In this invescigation a function of the form 

R * f <W D L- P r> E «- 1 

was sought, in which R * daily growth rate during the phenological interval, 

T * daily maximum temperature in °C, T = daily minimum temperature in °C, 
x m 

D ■ daily daylength in hours, and P * daily precipitation in millimeters. 

L IT 

Functional forms were selected for investigation, and parameters for these 
functions were estimated using generalized least squares techniques, to be 
discussed subsequently. For each phenological Interval Independent functions 


were sought. 
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Table 3: Typical Report of Dally Maximum OF POOR QUAl ITY 

and Minimum Temperatures 
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In general, two classes of functional forms were Investigated. The 
first was a polynomial of the form 


R * H, + H.T + H.T + H,D T + H C P + H,T 
1 2x 3m 4 L 5r ox 


+ H.T 2 + H q D t 2 + H q P 2 + H, n T D t 
7m 8 L 9r 10 xL 


+ H..T D t + H..P D. +R n TT + H..T P 
11 m L 12 r L 13 x m 14 x r 


+ H- -T P 
15 m r. 


The second was an exponential of the form 


R = exp(H + H.T + H.T + H.D T + H C P 
1 2x 3m 4 L 5r 

+.H.T 2 + H,T 2 + H Q D 2 + H 0 P 2 + H. rt T D T 
6x 7m 8 L 9r 10 x L 


+ H -T D. + H. 0 P D t + H n T T 
11 m L 12 r L 13 x m 


Eq. 2 


Eq. 3 


+ H. ,T P^ + H.-T P ). 
14 x r 15 m r 


A third functional form of interest is the "triquadratic" advocated by 
Robertson (1) for spring wheat, of the form 


R 


C a l (D L - V + a 2 (D L * a o> 2 ] ' D>l (T x ‘ V 
+ b 2 (T x ' V 2 + V T m - V + V T m - b o )2 l 


Eq. 4 


It should be noted, however, that this equation may be expanded, and a 
special case of the polynomial formulation results of the form 


R - »1 + H 2 T x + H 3 T ffl + H 4 D l + H 5 T x 2 + H 6 T ffl 2 

+ H 7 D L 2 + H 8 T x D L + V-A + "loV . 2 
+ h ii d l t » 2 + h i2 t x°l 2 + H n T A Z 


2 2 2 2 
+ H, , T D, + H. -T D t 
14 x L 15 m L 


Eq. 5 
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in which. 


H i ' VoVo ' Vo b 2 b o‘ + a l a O b 3 b O “ Vo b 4 b O 
‘ U 2 a o' b l b O + a 2 a 0 2b 2 b 0 2 ' a 2 a o' b 3 b O + a 2 a o‘ b * b O 
H 2 " 2a l a O b 2 b O _ a l“0 b l ‘ 2a 2 a O~ b 2 b O + a 2 a O b l 

H 3 * 2 V 0 b 4 b 0 ' *l a O b 3 " 2a 2 a o‘ b 4 b O + VoS 

H 4 * * ! b 2 b 0 2 ' a l b l b O + a l b 4 b O * a l b 3 b O + a l b 3 b O + 2a 2°O b l b O 


2 . 


- 2a 2 a 0 b 2 b 0 ‘ + 2a 2 a 0 b 3 b 0 ” 2a 2 a O b 4 b O 
H 5 " a 2 a O b 2 " l l a O b 2 

H b " a 2 a O b A “ Vo b A 

H 7 * a 2 b 2 b 0 2 ‘ a 2 b l b O + a 2 b A b 0 2 “ d 2 b 3 b o' 
H 8 " a l b l " 2a l b 2 b O - 2a 2 a O b l + 4a 2*O b 2 b O 

H 9 " a l b 3 " 2a l b A b O “ 2a 2 a O b 3 + 4a 2 a O b A b O 

H 10 " a l b 2 ’ 2a 2 a O b 2 

H U " a l b A ” 2a 2 a O b A 


H 12 " a 2 b l ‘ 2a 2 b 2 b O 


H 13 " a 2 b 3 “ 2a 2 b A b O 


*14 " a 2 b 2 


H 15 " a 2 b A 


F.q. b 


•SSSi**' 


It la recognized that this formulation results In an overparameterized ease 
when compared with Equation 4 In that 15 parameters result from the trans- 
formation to the polynomial form ns compared to only « In the original 


triquadratic. This is not a serious problem, however, if sufficient data are 
gathered to insure ample redundancy for the least squares estimation of 
parameters. Further, Robertson (1) suggests that depending upon the pheno- 
logical interval, simplified forms of Equation A be utilized in which some 
higher order terms are omitted. 

Least Squares Techniques 

For this investigation parameters for the above functions were estimated 
using the general constrained minimum least squares approach as published 
by Mikhail (3). The method, called "simultaneous adjustment of observations 
and parameters," allows maximum flexibility for investigations of alternate 
functional forms and provides an extremely powerful tool for determining 
parameters of these functions and estimating variances of the parameters 
obtained. The following discussion provides a brief summary of the technique. 
For a more detailed discussion, the reader is referred to Mikhail (3). 

Given a vector of n original observations 




2 * 



Eq. 7 


and a set of condition equations of the form 


F . L X _ <p 
r,l'n,l*u,l' r,l 


Eq . 8 


in which there are r equations, L is the vector of n adjusted observations, 
X is the vector of u parameters to be estimated, and ^ is the null vector, 
then the following may be stated. 

Since, in general, these equations will be nonlinear, a Taylor's series 


approximation may be written as 
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In which. 


_A V_ _B A_ B F° 

r,n n, 1 r,u u, l r, 1 


A 


.n. 



3 K 

- * 3X 


L°.X° 


Fq. ^ 


F.q. 10 


are the Jacobian matrices of the functions with respect to the observations 
and parameters, respectively, and evaluated at observed values and parameter 
estimates. Further, 


F . . F(l°,X 0 ) 
r,l — - — 


Eq. 11 


is the evaluation of the functions with observed values and parameter 
estimates, V is the vector of residuals to be applied to the observed 
quantities, and A is the vector of corrections to be applied to the 
parameters . 

The solution proceeds using the least squares condition in which 


v* _p v_ 

l , n n , n n , l 


Fq. 12 


is to be minimised. P is the weight matrix of the original observed data. 
If a priori est invites of variance are known, then 



n , n n , n 


F.q. 13 


in which 0 is the a priori variance/covariance matrix of the observed 
quantities. If this var lance/covar lance matrix is unknown, then V represents 
some assumed weight matrix, and may he written as 
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P - Eq. U 

in which ^ o represents some assumed cofactor matrix of relative variances 

L 

and covariances. Since, In general, the A Jacobian matrix in Equation 9 will 
not be square, and hence, the linearized condition equations cannot be 
solved directly for the residuals and substituted into the least squares 
condition of Equation 12, the more general method of constrained minimum 
is used. In this method, a vector of r Lagranglan multipliers is defined, 

- nk 1 .k 2 ,...k r J t . Eq. 15 

Then, the following scalar is minimized rather than A. 

- V* P V - 2k£(A V + R A - F°) Eq. 16 

Minimising by taking partials of this scalar with respect to the observa- 
tions and parameters and augmenting with the original linearized condition 
equations results in the following total set of normal equations: 


T 

euj 
1 


" _V_ ' 


" jl' 

n , n n , r n , u 


n, l 


n,i 

JL _1_ -1. 


A 


F° 

r,n r,r r,u 


r.l 


r.l 

JL JL _4_ 


_A_ 


JL. 

_u,n u,r tt s U M 


. u, 1 _ 


_ u.l . 


Eq. 17 


This set may be solved by partitioning with the following results: 

_A _ iT 1 T_ 
u, 1 U , U 11,1 


Eq. 18 


15 




» (A £ l o aV* (-B A + F°) 


Eq. 19 


1 ’ V £ *L 

in which. 


N 

u,u 



A ) B 


Eq. 20 


Eq. 21 


T 

u.l 


B C (A Sjq A^'V 


Eq. 22 


These equations are solved in an iterative manner until convergence. 
After convergence, the reference variance may be computed by 


V C P V 

s*" » 

o r-u 


Eq. 23 


Variance propagation techniques may be used to show that (3), 



Eq. 24 


yields an estimate of the a posteriori variance/covariance nuitrix for the 
parameters. 

The reference variance s~ represents a measure of the "goodness" of the 
fit, and may be used in F tests to determine the statistical significance in 
adding or deleting terms from the functional form under investigation, or 
in comparing different functional forms. 

A special case of the above procedure, called "variation of parameters," 
is more familiar to those who have used regression techniques. For this case. 
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only one observed quantity appears in each of the condition equations; that 
is, the number of observations is equal to the number of condition equations. 
In regression, for example, one observed value of the dependent variable, y, 
occurs in each equation. For this case, it may be seen that the A matrix in 
Equation 9 becomes identity. Equation 9 could, then, be solved directly for 
the residuals, this result substituted directly into Equation 12, and this 
sum of squares minimized with respect to the parameters. 

A more direct approach, however, is to substitute the identity matrix 
for A in Equations 21 and 22. The results are 


= B t PB 

Eq. 25 

c 

c 

1 

1 

1 

JL = b'pf 0 

u,l 

Eq. 26 


The solution then proceeds as before. 

Note that, with the variation of parameters method, whichever solutions 
technique is used, the method still does not depend upon linearity of the 
condition equations with respect to the parameters. Thus, this special case 
can still accomodate nonlinear functions through the use of the Taylor's 
series linearization and iterative techniques. 

For the problem at hand the function of Equation 1, formulated in terras 
of Equation 8, would be for each data point 

F - R - f (T ,T ,D ,P ) - 0. Eq. 27 

x m l r 

in which r - 1, n - 5, and u would depend upon the number of parameters 
required for the functional form selected. For the case using the general 
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least squares procedure, 

-r^T - R T T D t P r 
5»1 L x m L rj 

- H. H, ... H r. 
u,l [_ 1 2 uj 

If it were desired to use the variation of parameters, the independent 

variables T , T , D, , P would be assumed as fixed, or constants, for each 
x m L r 

data point and 


Eq. 28 
Eq. 29 


L 


1,1 


1 - R 


Eq. 30 


would result. 


IMPLEMENTATION 

Efforts were made to formulate and execute a data reduction system 
subject to the constraints of the data forms available and consistent with 
the data requirements for application of the least squares techniques dis- 
cussed in the previous section. Computer programs were written for the 
polynomial and exponential function forms using both the variation of para- 
meters and the more general least squares techniques. Special techniques 
were required in treating the Robertson's model. 

Data Reduction and Preparation 

Although the objective of the research was to model dally growth rates, 

R, as functions of environmental variables, no data could be obtained for 
growth rate on a daily basis. As may be seen from Table II, phonological data 
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reported consisted at best of dates at which the crop reached each pheno- 

logical stage. With this data, only an average growth rate could be computed, 

* 

of the form of 



in which N » the number of days in the i-th phenological interval for the 
location/year. This average growth rate represents an average, not only over 
the time period for the interval, but also a spatial average for the crop 
over the geographic region representing the crop reporting district. 

Environmental data were available on a daily basis. However, it would 
be inconsistent to apply daily environmental data in the least squares 
modeling process when only averaged growth rates were available. Therefore, 
in order to supply environmental data consistent with the growth rates avail- 
able, the following pre-processing steps were performed. 

First, within each crop reporting district several meterological 
reporting stations were selected. Up to five stations were used. Tempera- 
ture and precipitation data were recorded for each station on a daily basis 
throughout the reported phenological stage interval from the climatological 
records. Table 5 represents a typical data reduction form used for this 
purpose. 

This information was keypunched onto cards and input into a data re- 
duction program. This program computes the average position of the crop 
reporting district (latitude, longitude, elevation) and the associated 
statistics, converts the meterological data from English to SI units, 

averages meterological values (T ,T ,P ) for the location/year, and computes 

x ro r 









the associated variances. A listing of this data reduction program with 
output from a typical run is shown in Appendix A. For all of the pheno- 
logical intervals, with the exception of the emerge-joint interval, the re- 
sults for each location/year would be a single data point for input into the 
least squares program, containing averaged growth rate, maximum temperature, 
minimum temperature, daylength, precipitation, and positional information for 
that location/year. 

Special considerations were necessary for the emerge-joint interval. 
Because this interval stretches over the long winter dormancy period, no 
single averaged growth rate value could be considered to reasonably repre- 
sent the entire interval. For the growth rates the interval was broken into 
three sections, as depicted in Figure II. A fall emerge-dormant period was 
defined from the reported emergence date until average temperatures, by 
visual inspection of meterological records, appeared to fall below 4.5°C. 
(40°F.). This value of 40°F. was relatively arbitrary based to some extent 
upon informal conversations with agronomists at JSC and the USDA Research 
Laboratory at Fort Collins, Colorado (4,5). The crop was assumed to have 
completed one half of its interval growth in this period, and a growth rate 
of 


*e 



Eq. 32 


was assigned for this period. Over the winter dormancy period, Np, assumed 
to last until the average daily temperature again rose above 4.5°C. (40°F.) 


In the spring, a growth rate of zero was assumed. Lastly, a spring greenup- 
joint interval was defined, and a second average growth rate computed as 
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Eq. 33 


These rates, and the environmental data averaged on a monthly basis were 
used to define several data points for each location year in this interval. 

Computer Runs 

The functional forms of Equations 2 and 3 were programmed for parameter 

estimation in two modes. In the first, the method of variation of parameters 

was programmed. In using this method, . t must be assumed that the values 

of independent variables (T ,T ,D. ,P ) are perfectly known, and the sum of 

x m L r 

squares of residuals for the dependent variable, rate, is minimized. This 

technique most closely resembles the more familiar regression techniques 

common in statistics. A listing of this program, with typical output is 

shown in Appendix B. For the investigation, functional forms were selected 

and computer runs made using selected terms from Equations 2 and 3. For the 

polynomial, runs were made using only the constant term, using single variable 

linear equations in D , T , and T , using linear combinations of these 

L m x 

variables with and without precipitation, and various selected combinations 
of linear with interaction terms and linear with squared terms. A more 
limited subset of these combinations was run for the exponential function 
using the variation of parameters approach. 

In the second mode, the method of combined adjustment of observations 
and parameters was used. This represents the most generalized method of 
least squares curve fitting, and is henceforth referred to simply as 
"generalized least squares". Combinations of terms were run similar to those 
described above using this program. A listing of the program, together with 
output from a typical run, are shown in Appendix C. 
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Specific features of the programs which may be of interest include the 
fact that any functional form may be accomodated by changing a single statement 
in the EF function subprogram. Linearization is accomplished by numerical 
evaluation of partial derivatives using forward differences. For the gen- 
eralized least squares program (Appendix C) , even greater flexibility is pro- 
duced by the introduction of a parameter selection index, IPAR, which allows 
the inclusion or deletion of parameters within the functional form being 
investigated. 

The use of such numerical techniques, which allows great flexibility, is 
not without danger. It was found in some instances that convergence of the 
solution did not occur when it was to be expected due to the approximations 
inherent in these numerical techniques. These problems occurred only for the 
highly nonlinear functions and is not felt to be a serious problem for reasons 
to be discussed subsequently. 

The Robertson's Model 

As discussed previously, the "triquadratic" model advocated by Robertson 
has unich to recommend it. The model, Equation 4, can accommodate both non- 
linear effects and interactions with a minimum number of parameters. However, 
the model requires special considerations in estimating parameters. As 
pointed out by Robertson in his paper (1), the parameters for the model are 
not independent. An attempt to simultaneously estimate all of these para- 
meters using the generalized least squares technique verified this, for any 
attempt at such modeling results in a near-singular normal equation matrix. 

One solution to the problem, expansion of the equation into polynomial 
form, shown in Equation 5, is less than desirable. This expansion results in 
an overparameterized situation and the model loses its efficient character- 
istics. Therefore it is desirable to seek other methods for fitting an 
equation of the "trlquadratlc" form. 


IS 
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Robertson utilized a parameter transformation technique. By expansion 
of each multiplicative factor separately, polynomials result which have the 
same number of coefficients as original parameters within each factor. 
Estimation may then be carried out by holding the coefficients of the trans- 
formed daylength parameters, for example, and using regression to estimate 
transformed temperature parameters. These temperature parameters may then be 
held, and transformed daylength parameters found by regression. This 
successive substitution technique is cycled until no change in the parameters 
result, and the parameters are transformed back for the original equation 
form. The important point here is that the temperature parameters are 
estimated independently from those for daylength. 

The same technique may be implemented more simply using the generalized 
least squares techniques discussed previously. The procedure is identical, 
except no transformation is required, since the method is not limited to 
regression. The method proceeds as ioIIows. First, parameter approximations 
are assumed for the daylength factor. There are held as constants, and the 
general least squares algorithm are applied in an iterative manner until 
convergence for the temperature parameters. These are then held, and 
another computer run made to estimate new values for the daylength para- 
meters, using these generalized least squares algorithms. The process is 
repeated in successive computer runs until no change is noted in the 
parameters. 


TESTING PROCEDURES 

Two types of testing procedures are necessary to fully assess the 
adequacy of the models, to discriminate between equation forms, .u between 
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various parametric choices within each equation form. In the first, 
a -posteriori variances resulting from the least squares techniques may be 
compared using statistical tests. A more appropriate technique, however, is 
to apply the functions generated by least squares estimation co independent 
data, not included for parameter estimation. The first technique was fully 
applied during this investigation. The second method, which requires 
extensive data preparation, was applied over limited data sets only. 

Variance Comparisons 

Tables 6-9 depict the variances obtained a- posteriori , about the least 
squares fit for both the polynomial and exponential functional forms, obtained 
using the variation of parameters and the generalized least squares technique. 
The column headed "Terms Entered in Equation" refers to Equations 2 and 3, 
respectively. The degrees of freedom (d.f.) for each case represents the 
total number of parameters. selected for each fit. The variance values given, 
then, represent the sum of squares of the residuals for growth rates (R - R) 
divided by these degrees of freedom. 

In addition to the polynomial and exponential functions, the special 
case of the Robertson’s triquadratic model. Equation 4, was investigated. In 
Table 10, the variance about this equation are shown. 

Inspection of these tables reveals that, in general, the variances 
obtained using the generalized least squares technique exceed those obtained 
using the variation of parameters method. This is to be expected, since the 
variation of parameters method minimizes residuals only of growth rates, while 
the generalized least squares simultaneously minimize residuals for all 
observed quantities. Nevertheless, the parameters obtained from generalized 
least squares are considered the more reliable, sinca this technique more 




Stage 0-1 

Stage 1-2 

Stage 2-3 

Stage 3-4 

Stage 4-5 

Function 

Terms Entered 
in Equation* 

d.f. 

Var-. 

(x 10" 3 ) 

d.f. 

Var. 

(x 10" A ) 

d.f. 

Var . 

(x 10" 5 ) 

d.f. 

Var. 

(x 10~ 4 ) 

d.f. 

Var. 

(x 10" 4 ) 

Const. Only 

1 

30 

.6107 

2 38 

.3717 

43 

.6865 

38 

.8011 

43 

.3496 

Lin. Only 

1.4 

30 

.6107 

237 

.2795 

42 

.5147 

37 

.5516 

42 

.3576 

Lin. T Only 
m • 

1.3 

29 

.6099 

237 

.2447 

42 

.6560 

37 

.8173 

42 

.3510 

Lin. T x Only 

1.2 

29 

.5053 

237 

.2501 

42 

.6980 

37 

.8212 

42 

.3041 

T , D. , & T D. 
x* L x L 

1,2,4,10 

29 

.5053 

235 

.2260 

40 

. 3470 

35 

.5571 

40 

.3012 

T , T , & T T 
x m x m 

1-3,13 

27 

.5115 

235 

.2400 

40 

.6845 

35 

.8525 

40 

.3081 

Lin. T , T , D_ 
X m 9 T 

1-4 

28 

.5107 

235 

.2177 

40 

.4974 

35 

.5599 

40 

.3097 

Lin. T , T , D T , P 

1-5 

27 

.5296 

234 

.2186 

39 

.4960 

34 

.5583 

39 

.3070 

A 111 Lt L 

Lin. & Int. w/o ? r 

1-4,10,11,13 

27 

.5115 

232 

.2091 

37 

.2840 

32 

.5751 

37 

.3154 

Lin. & Sq. w/o ? r 

1-4, 6-8 

26 

.4389 

232 

.2077 

37 

.4782 

32 

.5566 

37 

.3062 

Lin. & Sq. w/ P^ 

1-9 

24 

.4683 

230 

• .2086 

35 

.4889 

30 

.5500 

35 

.3127 

All w/o P f 

1-4 j,6 3,10,11 ,13 

25 

. 3931 

229 

.1863 

34 

.2189 

29 

.5208 

34 

.2784 

Lin. & Int. w/ P 

r 

1-5.10-15 

24 

.4929 

228 

.2036 

33 

.1413 

28 

.5847 

33 

.3112 

All 

1-15 

21 

.4294 

224 

.1853 

29 

.1597 

24 

.5834 

29 

.2997 


*For Stage 0-1, all terms containing D ? were omitted 



Table 6: Variances About Polynomial 

Fit Using Variation of Parameters 




Stage 0-1 

Stage 1-2 

Stage 2-3 

Stage3- 4 

Stage 4-5 


Function 

Terms Entered 
in Equation* 

j f var.- 

d * ' (x 10 3 ) 

Var . , 

d,f * (x 10 -4 ) 

j c r ■ , 

d,f ’ (x 10" 5 ) 

a f Var.. 

d,f * (x 10" 4 ) 

, r Var., 

d,f * (x 10" 4 ) 




Const. Only 

1 

30 

.6107 

238 

.3718 

43 

.6865 

38 

.8011 

43 

.3496 

Lin. T , T , D. 
x m L 

1-4 

28 

.4951 

235 

.2152 

40 

.4342 

35 

.5460 

40 

.3116 

Lin. T , T , D_ , P 
x m L - r 

1-5 

27 

.5133 

234 

.2161 

39 

.4297 

34 

.5428 

39 

.3064 

Lin. h Int. w/o P f 

1-4,10,11,13 

27 

.4965 

CM 

cn 

CM 

.1865 

37 

.1546 

32 

.5560 

37 

.3210 

Lin. & Sq. w/o P r 

1-4, 6-8 

26 

.4329 

232 

.1631 

37 

.2985 

32 

.5764 

37 

.3036 

Lin. & Sq. w/ P 

r 

1-9 

24 

.4631 

230 

.1640 

35 

.2910 

30 

.5390 

35 

.3094 

All w/o P f 

1-4,6-8,10,11,13 

25 

.3603 

229 

.1285 

34 

. 1351 

29 

.5012 

34 

.2723 

Lin. A Int. w/ P 

r 

1-5,10-15 

24 

.4651 

228 

.1534 

33 

.1034 

28 

.5711 

33 

.3239 

All 

1-15 

21 

. 3856 

224 

. 1252 

29 

.1108 

24 

.5710 

29 

.2945 


*For Stage 0-1, all terms containing were omitted 


Table 7: Variances About Exponential 

Fit Using Variation of Parameters 





Stage 0-1 

Stage 1-2 

Stage 2-3 

St ape 3-4 

Stage 4-5 

Function 

Terms Entered , f Var io 

in Equation* (x 10 ) 

, c Var., 

d-f - (x 10 A ) 

j f Var.,. 

d,f * (x 10" 5 ) 

a c Var. , 

d,f * (x 10" 4 ) 

a t Var., 

d,f# (x 10 _A ) 



Const. Only 

1 

30 

.6107 

238 

.3717 

43 

.6865 

38 

.8011 

43 

.3496 

Lin. D Only 

1,4 

30 

.6107 

237 

.3638 

42 

.4925 

37 

.5516 

42 

.3576 

Lin. T Only 
m 

1,3 

29 

.6099 

237 

.3104 

42 

.6560 

37 

.8173 

42 

.3510 

Lin. T Only 

X 

1,2 

29 

.5053 

237 

.3187 

42 

.6980 

37 

.8212 

42 

.3041 

T , D_ , & T D. 
x L x L 

1,2,4,10 

29 

.5053 

235 

.3058 

40 

.5273 

35 

.5574 

40 

.3029 

Lin. T , T , D_ 
x m L 

1-4 

28 

.5107 

235 

.3035 

40 

.4549 

35 

.5599 

40 

.3079 

Lin. T , T , D , P 
x m* L r 

1-5 

27 

.5296 

234 

.3042 

39 

.4689 

34 

.5583 

39 

.3070 

Lin. & Int. w/o P^ 

1-4,10,11,13 

27 

.5159 

232 

.3050 

37 

.3521 

32 

.5837 

37 

.3233 

Lin. & Sq. w/o P^ 

1-4, 6-8 

26 

.4429 

232 

.2663 

37 

.5173 

32 

.5568 

37 

.3179 

Lin. Sq. w/ P 

r 

1-9 

24 

.4797 

230 

.2668 

35 

.5259 

30 

.5712 

35 

.3150 

All w/o P 

r 

1-4,6-8,10,11,13 

25 

.4183 

229 

.2522 

34 

.3025 

29 

.5262 

34 

.3292 

Lin. & Int. w/ P 

r 

1-5,10-15 

24 

.5023 

228 

.3459 

33 

.2545 

28 

(.6049) 

33 

.3488 

All 

1-15 

21 

.5465 

224 

.2611 

29 

.2888 

24 

(.8071) 

29 

.3397 


*For Stage 0-1, all terms containing were omitted 
( ) Solution did not converge 


Table 8: Variances About Polynomial 

Fit Using Generalized Least Squares 




Stage 0-1 

Stage 1-2 

Stage 2-3 

Stage 3-4 

Stage 4-5 

Function 

Terms Entered 
in Equation* 

d.f. 

Var. ^ 
(x 10" J ) 

d.f. 

Var. , 
(x 10"*) 

d.f. 

Var. - 

(x io~h 

d.f. 

Var. . 
(x 10 4 ) 

d.f. 

Var . , * 
(x 10 *) 


Const. Only 

1 

30 

.6107 

238 

.3717 

43 

.6865 

38 

.8011 

43 

.3496 

Lin. Only 

1,4 

30 

.6107 

237 

.3824 

42 

.4923 

37 

.5436 

42 

.3755 

Lin. T Only 

X 

1,2 

29 

.4906 

237 

.3345 

42 

.6976 

37 

.8211 

42 

.3054 

T , D , & T D_ 
x* L x L 

1,2,4,10 

29 

.4906 

235 

.2976 

40 

.5599 

35 

.5450 

40 

3050 

Lin. T , T , D, 
x* m* L 

1-4 

28 

.4965 

235 

.3123 

40 

.4350 

35 

.5461 

40 

.3119 

Lin. T , T , D , P 
x* nr L* r 

1-5 

27 

.5148 

234 

.3118 

39 

.4263 

34 

.5437 

39 

.3074 

Lin. & Int. w/o P 

r 

1-4,10,11,13 

27 

.5040 

232 

.3610 

37 

(.6389) 

32 

.5625 

37 

.3358 

Lin. & Sq. w/o P^ 

1-4, 6-8 

26 

.4494 

232 

.2317 

37 

.4588 

32 

.5561 

37 

.3060 

Lin. & Sq. w/ P^ 

1-9 

24 

.5274 

230 

.23;6 

35 

.4345 

30 

.5702 

35 

.3128 

All w/o P r 

1-4,6-8,10,11,13 

25 

.3781 

229 

.2515 

34 

.2549 

29 

.5116 

34 

.2887 

Lin. & Int. w/ P 

r 

1-5,10-15 

24 

.4843 

228 

.3016 

33 

.1562 

28 

.5965 

33 

(.4065) 

All 

1-15 

21 

.4501 

224 

.2418 

29 

(.2565) 

24 

(.8071) 

29 

(.3467) 


*For Stage 0-1, all terms containing D^were omitted 
( ) Solution did not converge 


Table 9: Variances About Exponential 

Fit Using Generalized Least Squares 
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clearly recognizes the physical situation of randomness in the independent 

variables T , T , D T , P . 

x m L r 

Table 10: Variances About "Triquadratic" Fit 

Using Generalized Least Squares 


Parameters Degrees of 

Stage Estimated Freedom Variance 


0-1 

b 0* b l ,b 2' b 3* b 4 

26 

.4429 x 10 -3 




-4 

1-2 

a 0' a l ' b 0 ,b l ,b 2 ,b 3 ,b 4 

232 

.2987 x 10 




-3 

2-3 

a 0 ,a l* b 0* b l' b 3 

39 

.4465 x 10 

3-4 

a 0' a i ,b 0 ,b l* b 3 

34 

,5794 x 10~ 4 




-3 

4-5 

a 0’ a !' b 0 ,b l* b 3 

39 

.3192 x 10 


These variances may be compared by forming ratios, and testing with the 
statistical F-test. This was done for a variety of combinations, in an 
attempt to assess the following: 

1. ) in general, is there any advantage in using the exponential 

function over the more easily applied polynomial; 

2. ) does the inclusion of precipitation make any statistically 

significant difference in the model; 

3. ) is the Robertson's "triquadratic" model with it advantages, a 

viable alterantive function to the polynomial and exponential 

equations for winter wheat; 

4. ) within each equation form, which terms are most significant? 

In Tables 11 and 12, ratios of variance from the polynomial over those 
for the exponential are shown, for comparable cases. In each instance, the 
number of parameters in each model, and hence the degrees of freedom, are 







Stage 0-1 

Stage 1-2 

Stage 2-3 

Stage 3-4 

Stage 4-5 

Function 

Terms Entered . . Var. 

in Equation** * * (Ratio) 

j c Var. 

* * (Ratio) 

j f Var. 

* ' (Ratio) 

j f Var. 

* • (Ratio) 

o f Var. 

a,r * (Ratio) 



Const. Only 

1 

30 

1.000 

238 

1.000 

43 

1.000 

38 

1.000 

43 

1.000 

Lin. T , T , D. 
x m* L 

1-4 

28 

1.031 

235 

1.012 

40 

1.146 

35 

1.025 

40 

0.994 

Lin. T , T , D, , P 
x m* L r 

1-5 

27 

1.032 

234 

1.012 

39 

1.154 

34 

1.028 

39 

1.002 

Lin. 4 Int. w/o P 

r 

1 

«0 

o 

V 

u> 

27 

1.C30 

232 

1.121* 

37 

1.837* 

32 

1.034 

37 

0.983 

Lin. 4 Sq. w/o P^ 

1-4, 6-8 

26 

1.014 

232 

1.273* 

37 

1.602* 

32 

0.966 

37 

1.009 

Lin. 4 Sq. w/ P^ 

1-9 

24 

1.011 

230 

1.272* 

35 

1.680* 

30 

1.020 

35 

1.011 

All w/o P r 

1-4,6-8,10,11,13 

25 

1.091 

229 

1.450* 

34. 

", 1.620* 

29 

1.039 

34 

1.022 

Lin. 4 Int. w/ P 

r 

1-5,10-15 

24 

1.060 

228 

1.327* 

33 

1.366* 

28 

1.024 

33 

0.961 

All 

1-15 

21 

1.114 

224 

1.480* 

29 

1.441* 

24 

1.022 

29 

1.018 


^Indicates significance at a - 0.25 
**For Stage 0-1, all terms containing were omitted 


Table 11: Patios of Variance About 

Fit (Polynomial/Exponential) Using Variation of Parameters 


Stage 0-1 
Rbt . df=26 

Stage 1-2 
Rbt. df=232 

Stage 2-3 
Rbt. df-39 

Stage 3-4 
Rbt. d f» 34 

Stage 4-5 
Rbt. df-39 

„ , . . _ . Terms Entered n . 

Polynomial Function . ^ . d.f. Ratio 

* in Equation 

d.f. Ratio 

d.f. Ratio 

d.f. Ratio 

d.f. Ratio 



Lin. T , T , D T 

1-4 

28 

0.867 

235 

0.984 

40 

0.982 

35 

1.035 

40 

1.031 

x ra* L 
Lin. & Int. w/o 

1-4,10,11,13 

27 

0.858 

232 

0.979 

37 

1.286** 

32 

0.993 

37 

0.987 

Lin. & Sq. w/o P f 

1-4, 6-8 

26 

1.000 

232 

1.122** 

37 

0.863 

32 

1.040 

37 

1.004 

All w/o P f 

1-4,6-8,10,11,13 

25 

1.059 

229 

1.184** 

34 

1.476** 

29 

1.101 

34 

0.970 


Exponential Function 


Lin T , T ,D t 
x m L 

1-4 

28 

0.829 

235 

0.956 

40 

1.026 

35 

1.061 

40 

1.023 

Lin. & Int. w/o P 

r 

1-4,10,11,13 

27 

0.876 

232 

0.827* 

37 

(0.699) 

32 

1.030 

37 

0.951 

Lin. & Sq. w/o P^ 

1-4, 6-8 

26 

3.986 

232 

1.289** 

37 

0.973 

32 

1.042 

37 

1.043 

All w/o P f 

1-4,6-8,10,11,13 

25 

1. 171 

229 

1.187** 

34 

1.752** 

29 

1.132 

34 

1.106 


♦Indicates significance in Robertson over other at a b 0.25 
♦♦Indicates significance in other over Robertson at a * 0.25 


Table 13: Variance Ratios for Comparison 

of Robertson's Model With Other Functions 


Function 


Terms Entered 
in Equation** 


Stage 0-1 Stage 1-2 Stage 2-3 Stage 3-4 Stage 4-5 


x 10 ') (x 10 (x 10^) 


Const. Only 

1 

30 


238 

1.000 

43 

1.000 

38 

1.000 

43 

1.000 

Lin. Only 

1.4 

30 


237 

0.951 

42 

1.000 

37 

1.015 

42 

0.952 

Lin. T^ Only 

1,2 

29 

1.030 

237 

0.952 

42 

1.001 

37 

1.000 

42 

0.996 

V V 4 T * D L 

1,2,4,10 

29 


235 

1.028 

40 

0.942 

35 

1.023 

40 

0.993 

Lln - V V D L 

1-4 

28 

1.029 

235 

0.972 

40 

1.046 

35 

1.025 

40 

0.987 

Lin. T . T , D_, P_ 
x nr L r 

1-5 

27 

1.029 

234 

0.976 

39 

1.100 

34 

1.027 

39 

0.999 

Lin. & Int. w/o P f 

1-4,10,11,13 

27 

1.024 

232 

0.845 

37 

(0.551) 

32 

1.038 

37 

0.963 

Lin. & Sq. w/o P f 

1-4, 6-8 

26 

0.986 

232 

1.149* 

37 

1.128 

32 

1.001 

37 

1.039 

Lln. & Sq. w/ P 

r 

1-9 

24 

0.910 

230 

1.152* 

35 

1.210 

30 

1.002 

35 

1.007 

All w/o P 

r 

1-4,6-8,10,11,13 

25 

1.106 

229 

1.003 

34 

1.187 

29 

1.029 

34 

1.140 

Lln. & Int. w/ P 

1-5,10-15 

24 


228 

1.147* 

33 

1.629* 

28 

(1.014) 

33 

(0.858) 

r 

All 

1-15 

21 

1.214 

224 

1.080 

29 

(1.126) 

24 

(1.000) 

29 

(0.980) 


*Indicates significance at a * 0.25 
**For Stage 0-1, all terms containing were omitted 
( ) Solution did not converge 


Table 12: Ratios of Variance About 

Fit (Polynomial/Exponential) Using Generalized Least Squares 
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identical for the numerator and denominator. 

Inspection of these tables reveals that the fit using exponential 
functions was significantly better than that using polynomials in only a few 
stages, and only for the higher-order cases. As will be discussed in the 
next section, the application of higher-order polynomials or exponentials to 
daily test data is dangerous. In general, then, it may be assumed that the 
use of exponential functions will not generally improve the least squares 
fit. 

The next assessment to be made concerns the inclusion of precipitation 

in the model. Visual inspection of Tables 6-9 reveals that the addition of 

precipitation terms over any of the existing models containing T^, T^, and 

D , in no instance significantly bettered the variance of the fit. It may 
L 

then be ten3tively stated that the inclusion of precipitation within a 
polynomial or exponential may be expected to have little effect (or a 
detrimental effect) upon the least squares fit. 

In evaluating Robertson's "triquadratic" model, a similar procedure, 
utilizing variance ratios may be used. The model is efficient in that non- 
linear and interaction effects may be accommodated with a minimum of 
parameters, and the computations required in later application of the 

equation are relatively simple and fast. It was found that, for the data 

2 

set used, no convergence of the solution could be obtained if the term 
was included. Thus, the cases run Included only linear terms. For the 

plant-emerge interval the entire daylength factor was set equal to one, 
recognizing that daylength could have no effect during that interval. The 
temperature factor Included the terms recommended by Robertson in his 
investigations of spring wheat (1). 
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Table 13 represents the ratios of variances obtained from the 
Robertson's fit compared with those from the appropriate cases from the 
polynomial and exponential functions generated using generalized least 
squares. In all cases, the variance from the "triquadratic" is the numerator. 
Few conclusions may be drown from these comparisons. In general, Robertson's 
model appears significantly better in some instances, and othir equation 
forms in others. One comment should be made, however. In the stages where 
significance appears, the parameter values converged to in the Robertson's 
fit did not appear "logical", in terms of the known phcnological behavior 
of winter wheat. The model is highly nonlinear and hence is dependent to a 
great extent on good parameter estimates to start the iterative runs. It is 
this investigator's opinion that the Robertson's model continues to show 
considerable promise, and work is continuing on this technique. For the 
present, however, the use of polynomial or exponential function forms is 
recommended for application and further testing. 

The last assessment performed using variances from the fits was in- 
tended to throw light on which terms were most significant within each 
equation form. Tables 14-18 depict variance ratios obtained from the various 
polynomial configurations considered, using the variation of parameters 
technique, for each stage. Tables 19-23 show similar information for the 
exponential functions fit using variation of parameters. Tables 24-28, and 
29-33 depict similar information for the polynomial and exponential func- 
tions, respectively, using the generalized least squares technique. 

The first case run, that of "Constant Only" represents the simplest 
possible assumption - that winter wheat growth rates are unaffected by 
environmental variations, and the crop matures at a constant rate within 
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in Equation 
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1 

1 
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1 

H 

rH 

rH 

rH 

H 

H 

rH 

rH 

rH 

H 

*H 

rH 


d.f. 30 30 29 29 29 28 27 27 27 26 25 24 24 21 


Const. Only 

30 

1.000 













Lin. , Only 

30 

1.000 

1.000 












Lin. , T Only 
in 

29 

1.001 

1.001 

1.000 











Lin. , T Only 

29 

1.208 

1.208 

1.207 

1.000 










V D L ‘ T x D L 

29 

1.208 

1.208 

1.207 

1.000 

1.000 









Lin., T x> T m , D l 

28 

1.194 

1.194 

1.194 

0.989 

0.989 

1.000 








T T , T T 
x m x m 

27 

1.196 

1.196 

1.192 

0.988 

0.988 

0.998 

1.000 







Lin., T x , T m , D l , P r 

27 

1.153 • 

1.153 

1.152 

0.954 

0.954 

0.964 

0.966 

1.000 






Lin. 6 Int. w/o P r 

27 

1.194 

1.194 

1.192 

0.988 

0.988 

0.999 

1.000 

1.035 

1.000 





Lin. & Sq. w/o 

26 

1.391* 

1.391 

1.389 

1.151 

1.151 

1.164 

1.165 

1.207 

1.165 

1.000 




All w/o P 

r 

25 

1.554 

1.554 

1.551 

1.285 

1.285 

1.299 

1.301 

1.347 

1.301 

1.117 

1.000 



Lin. & Sq. w/ P r 

24 

1.304 

1.304 

1.302 

1.079 

1.079 

1.091 

1.092 

1.131 

1.092 

0.937 

0.839 

1.000 


Lin. & Int. v/ P 

r 

24 

1.239 

1.239 

1.237 

1.025 

1.025 

1.036 

1.038 

1.074 

1.038 

0.890 

0.798 

0.950 

1.000 

All 

21 

1.422 

1.422 

1.420 

1.177 

1.177 

1.189 

1.191 

1.233 

1.191 

1.022 

0.915 

1.091 

1.148 1.000 


♦Indicates significance at a 3 0.25 


Table 14: Ratios of Variance About Polynomial Fit 

Using Variation of Parameters, Stage 0-1 


Terms Entered 
in Equation 








o 

rH 

0k 

H 



0 

»H 

H 

0k 

o 

•H 

00 

1 

NO 


o 

i-i m 

• ^H 

00 1 

1 o 

NO «~l 





m 


CJ 

m 


iTt 

'S' 

'S 

O' 

'S m «— i 



rH 

H 

rH 

tH 

H 

H 

*H 

fH 

fH 

rH 

•H 

1 1 1 

rH rH H 

Funct ion 

d.f. 

2 38 

237 

237 

237 

235 

235 

235 

2 34 

232 

232 

230 

229 228 224 

Const. Only 

238 

1.000 












Lin. Only 

237 

1 . 330 * 

1.000 











Lin. T Only 
m 

237 

1.519 

1 . 142 * 

1.000 










Lin. T x Only 

237 

1.486 

1.118 

0.978 

1.000 









T.Dl.TD, 
x L x L 

235 

1.645 

1.237 

1.083 

1.107 

1.000 








T , T , T T 
k m x m 

235 

1.549 

1.165 

1.020 

1.042 

0.942 

1.000 







Lin. T , T , D 
x* m* L 

235 

1.707 

. 1.284 

1 . 124 * 

1.149 

1.038 

1.102 

1.000 






Lin. T , T , D , P 
x m l r 

234 

1.700 

1.279 

1.119 

1.144 

1.034 

1.098 

0.996 

1.000 





Lin. & Int. w/o P 

r 

232 

1.778 

1.337 

1.170 

1.196 

1.081 

1.148 

1.041 

1.045 

1.000 




Lin. 6 Sq. w/o P f 

232 

1.790 

1.346 

1.178 

1.204 

1.088 

1.156 

1.048 

1.052 

1.007 

1.000 



Lin. (m Sq. w/ P f 

2 30 

1.782 

1.340 

1.173 

1.199 

1.083 

1.151 

1.044 

1.048 

1.002 

C . 996 

1.000 


All *w/o P 

r 

229 

1.995 

1.500 

1.313 

1.342 

1.213 

1.288 

1 . 169 * 

1.173 

1.122 

1.115 

1.120 

1.000 

Lin. 6 Int. w / P 

r 

228 

1.826 

1.373 

1.202 

1.228 

1.110 

1.179 

1.069 

1.074 

1.027 

1.020 

1.025 

0.915 1.000 

All 

224 

2.006 

1 . 508 

1.321 

1. 350 

1.220 

1.295 

1.175 

1.180 

1.128 

1.121 

1.126 

1.005 1.099 1.000 


*Tndicates significance at a = 0.25 


Table 15: Ratios of Variance About Polynomial 

Fit Using Variation of Parameters, Stage 1-2 


Terms Entered 
in Equation 


m 







* 



c 







rH 



rH 

m 


o 




rH 




rH 


rH 





oo 


00 

i 






c 

i 


1 

o 



i-H 



»H 

o 


vC 

rH 



* 



* 

+ 


* 

* 

m 

CM 

m 


m 

<r 

<T 

O' 

<r 

m 

H 

H 

H 

rH 

rH 

rH 

rH 

rH 

rH 

H 

rH 



Function 

d. f. 

43 

42 

42 

42 

40 

40 

40 

39 

37 

37 

35 

34 

33 

29 


Const. Only 

43 

1.000 














Lin . 

Only 

42 

1.334* 

1.000 













Lin. 

T Only 
m 

42 

1.046 

0.785 

1.000 












Lin. 

T Only 
x 7 

42 

0.984 

0.737 

0.940 

1.000 











V 

D t & T D, 
L x L 

40 

1.978 

1.483* 

1.890 

2.012 

1.000 










V 

T & T T 
m x m 

40 

1.003 

0.752 

0.958 

1.020 

0.507 

1 . 000 









Lin. 

T , T , D 
x m L 

40 

1.380 

1.035 

1.319 

1.403 

0.698 

1.376 

1.000 








Lin. 

T , T , D. , P 
x* m* L* r 

39 

1.384 

1.038 

1.323 

1.403 

0.700 

1.380 

1.003 

1.000 







Lin . 

& Int. w/o P 

r 

37 

2.417 

1.812 

2.310 

2.458 

1.222 

2.410 

1.751 

1.746 

1.000 






Lin. 

& Sq. w/o P 

r 

37 

1.436 

1.076 

1.372 

1.460 

0.726 

1.431 

1.040 

1.037 

0.594 

1.000 





Lin. 

& Sq. w/ P 

r 

35 

1.404 

1.053 

1.342 

1.428 

0.710 

1.400 

1.017 

1.015 

0.581 

0.978 

1.000 




All 

w/o P 

r 

34 

3.136 

2.351 

2.997 

3. 189 

1.585* 

3.127 

2.272 

1.266 

1.297 

2.185 

2.225 

1.000 



Lin . 

& Int. w/ P 

r 

33 

4.858 

3.643 

4.643 

4.940 

2.456 

4.844 

3.520 

3.510 

2.010 

3.384 

3.460 

1.549* 

1.000 


All 


29 

4.299 

3.223 

4.108 

4.371 

2.173 

4.286 

3. 115 

3.106 

1.778 

2.004 

3.061 

1.371 

0.885 

1.000 


*Indicates significance at a = 0.25 


Table 16: Ratios of Variance About Polynomial 

Fit Using Variation of Parameters, Stage 2-3 


Terms Entered 
in Equation 


m 

H 

I 








rH 

<* 

CO 

rH 



o 

rH 

,6-8 


,6-8 

o 

rH 





CO 

CN 

CO 

co 

i 


m 

i 

<r 

<r 

i 

1 

| 

m 

1 



rH 

rH 

rH 

r-H 

rH 

i 

rH 

r~H 

i 

rH 

r-H 

i 

rH 

rH 

rH 

rH 


Function 

d. f. 

38 

37 

37 

37 

35 

35 

35 

34 

32 

32 

30 

29 

28 


Con^L. Only 

38 

1.000 













Lin. Only 

37 

1.452* 

1.000 












Lin. T Only 
m J 

37 

0.980 

0.675 

1.000 











Lin. T^ Only 

37 

0.976 

0.672 

0.995 

1.000 










T , D , & T D 
x L x L 

35 

1.4 38 

1.000 

1.467 

1.474 

1.000 









T , T , & TT 
x m X m 

35 

0.940 

0.647 

0.959 

0.963 

0.653 

1.000 








Lin. T , T , D. 
x m L 

35 

1.431 

0.985 

1.460 

1.467 

0.995 

1.523 

1.000 







Lin. T T D , P 

x m L r 

34 

1.435 

0.988 

1.464 

1.471 

0.998 

1.527 

1.003 

1.000 






Lin. & Int. w/o P 

32 

1.393 

0.959 

1.421 

1.428 

0.969 

1-.482 

0.975 

0.971 

1.000 





Lin. & Sq . w/o P 

r 

32 

1.439 

0.991 

1.468 

1.475 

1.001 

1.532 

1.006 

1.003 

1.033 

1.000 




Lin. & Sq. w/ P 

r 

30 

1.457 

1.003 

1.486 

1.493 

1.013 

1.550 

1.018 

1.015 

1.046 

1.012 

1.000 



All w/o P 

r 

29 

1.538 

1.059 

1.569 

1.577 

1.070 

1.637 

1.075 

1.072 

1.104 

1.069 

1.056 

1.000 


Lin. & Int. w/ P 

r 

28 

1.370 

0.943 

1.398 

1.404 

0.953 

1.458 

0.958 

0.955 

0.984 

0.952 

0.941 

0.891 

1.000 

All 

24 

1.373 

0.945 

1.399 

1.408 

0.955 

1.461 

0.960 

0.957 

0.986 

0.954 

0.943 

0.893 

1.002 


LO 

f- ( 

I 


24 


* Indicates significance at a = 0.25 


Table 17: Ratios of Variance About Polynomial 

Fit Using Variation of Parameters, Stage 3-4 


ORIGINAL PAGE II* 
OF POOR QUALITY 


Terms Entered 
in Equation 



Function 

d.f. 

43 

42 

42 

42 

40 

40 

40 

39 

37 

37 

35 

34 

33 29 


Const. Only 

43 

1.000 













Lin. Only 

42 

0.978 

1.000 












Lin. T Only 
m 

42 

0.996 

1.019 

1.000 











Lin T Only 

X 

42 

1.150 

1.176 

1.154 

1.000 










T , D , & T D 
x L x L 

40 

1.161 

1.187 

1.165 

1.010 

1.000 









T , T , & T T 
x m x m 

40 

1.135 

1.161 

1.139 

0.987 

0.978 

1.000 








Lin. T , T , D. 
x in L 

40 

1.129 

1.155 

1.133 

0.982 

0.973 

0.995 

1.000 







Lin. T , T D P 

x m l r 

39 

1.139 

1.165 

1.14 3 

0.991 

0.981 

1.004 

1.009 

1.000 






Lin. 4 Int. w/o P 

r 

37 

1.108 

1.134 

1.113 

0.964 

0.955 

0.977 

0.982 

0.973 

1.000 





Lin. & Sq. w/o P 

r 

37 

1.142 

1.168 

1.146 

0.993 

0.984 

1.006 

1.011 

1.003 

1.030 

1.000 




Lin. & Sq. w/ P 

r 

35 

1.118 

1.144 

1.122 

0.972 

0.963 

0.985 

0.990 

0.982 

1.009 

0.979 

1.000 



All w/o P 

r 

34 

1.256* 

1.284 

1.261 

1.092 

1.082 

1.107 

1.112 

1.103 

1.133 

1.100 

1.123 

1.000 


Lin. L Int. w / P 

r 

33 

1.123 

1.149 

1.128 

0.977 

0.968 

0.990 

0.995 

0.986 

1.013 

0.984 

1.005 

0.895 

1.000 

All 

29 

1.166 

1.193 

1.171 

1.015 

1.005 

1.028 

1.03 3 

1.024 

1.052 

1.022 

1.043 

0.929 

1.038 1.000 


*Indicates significance at a = 0.25 


Table 18: Ratios of Variance About Polynomial 

Fit Using Variation of Parameters, Stage 4-5 


Terms Entered 
in Equation 







m 




ro 


i— t 




r— H 

ON 

1 






m 



r>. 



iH 


CO 

* 


1 


m 

rH 


VD 

in 

in 

m 

co 

m 


m 

m 

■ 




rH 

rH 

rH 

rH 

r —1 


— 





Function 

d.f. 

30 

28 

27 

27 

26 

25 

24 

24 

21 

Const. Only 

30 

1.000 









Li "' T x’ T m* °L 

28 

1.233 

1.000 








Lin - T x’ T m* D L> P r 

27 

1.190 

0.965 

1.000 







Lin. & Int. w/o P f 

27 

1.230 

0.997 

1.034 

1.000 






Lin. & Sq. w/o P^. 

26 

1.411* 

1.144 

1.186 

1.147 

1.000 





All w/o P f 

25 

1.695 

1.374 

1.425 

1.378 

1.201 

1.000 




Lin. & Sq. w/ P r 

24 

1.319 

1.069 

1.108 

1.072 

0.935 

0.778 

1.000 



Lin. & Int. w/ P r 

24 

1.313 

1.064 

1.104 

1.068 

0.931 

0.775 • 

0.996 

1.000 


All 

21 

1.584 

1.284 

1.331 

1.288 

1.223 

0.934 

1.201 

1.206 

1.000 


indicates significance at a = 0.25 


Table 19: Ratios of Variance About Exponential Fit 

Using Variation of Parameters, Stage 0-1 


-3,5-7,9,13-15 


Terms Entered 
in Equation 


m 


•a- 

I 


m 

I 


l 



Function d.f. 238 235 234 232 232 230 229 228 224 


Const. Only 

238 

1.000 


Lin. T , T , D 
x m L 

235 

1.727* 

1.000 

Lin. T , T , D , P 
x* m* L* r 

234 

1.720 

0.995 

Lin. & Int. w/o P f 

232 

1.993 

1.154* 

Lin. & Sq. w/o P^ 

232 

2.279 

1.319 

Lin. & Sq. w/ P f 

230 

2. 26r6 

1.312 

All w/o P^ 

229 

2.893 

1.675 

Lin. & Int. w/ P 

r 

228 

2.423 

1.403 

All 

224 

2.969 

1.718 


1.000 






1.159 

1.000 





1.325 

1.143* 

1.000 




1.318 

1.137 

0.995 

1.000 



1.682 

1.451 

1.269* 

1.276 

1.000 


1.409 

1.216 

1.063 

1.069 

0.838 

1.000 

1.726 

1.490 

1.303 

1.310 

1.026 

1.225 


1.000 


♦Indicates significance at a = 0.25 


ro 


Table 20: Ratios of Variance About Exponential 

Fit Using Variation of Parameters, Stage 1-2 


Terms Entered 
in Equation 


I 


ro 



«k 



O 



rH 



rH 



rH 



0k 

in 


» 

00 


00 

H 


O 

| 


1 

1 


H 

v£> 


NO 

rH 


r 

•» 


• 

* 

in 

| 

<r 

l 

ON 

| 

vr 

■ 

m 

1 

rH 

| 

rH 

rH 

1 

rH 

V 

f— 1 

1 

»H 

H 


Function d.f. 43 40 39 37 37 35 34 33 29 


Const. Only 

43 

1.000 


Lin. T , T , D 
x* m* L 

40 

1.581* 

1.000 

Lin. T , T , D , P 
x m L r 

39 

1.597 

1.010 

Lin. 4 Int. w/o P^ 

37 

4.441 

2.808* 

Lin. & Sq. w/o P f 

37 

2.300 

1.455 

Lin. & Sq. w/ P r 

35 

2.359 

1.492 

All w/o P f 

34 

5.081 

3.214 

Lin. & Int. w/ P r 

33 

6.639 

4.199 

All 

29 

6.196 

3.919 


1.000 






2.779 

• 1.000 





1.440 

0.518 

1.000 




1.477 

0.531 

1.026 

1.000 



3.181 

1.144 

2.209 

2.154 

1.000 


4.156 

1.495* 

2.887 

2.814 

1.306 

1.000 

3.878 

1.395 

2.694 

2.626 

1.219 

0.933 


1.000 


♦Indicates significance at a * 0.25 


O Q 
►*3 g 

g2 

i- 


-t* 

u> 




. 




Table 21: Ratios of Variances About Exponential 

Fit Using Variation of Parameters, Stage 2-3 


Terms Entered 
in Equation 


CO 


m 

I 


*3- 

I 


00 

I 

V0 


sf 

I 


cr> 

l 


co 


oo 

I 

VO 


I 


LO 


I 

o 


m 

tH 

I 


Function d.f. 38 35 34 32 32 30 29 28 24 


Const. Only 

38 

1.000 








Lin. T T , D 
x m l 

35 

1.467* 

1.000 







Lin T , T , D , P 
x* :o* L* r 

34 

1.476 

1.006 

1.000 





• 

Lin. & Int. w/o P f 

32 

1.441 

0.982 

0.976 

1.000 





Lin. & Sq. w/o P r 

32 

1.390 

0.947 

0.942 

0.965 

1.000 




Lin. & Sq. w/ P r 

30 

1.486 

1.103 

1.007 

1.032 

1.069 

1.000 



All w/o P^ 

29 

1.598 

1.089 

1.083 

1.109 

1.150 

1.075 

1.000 


Lin. & Int. w/ P f 

28 

1.403 

0.956 

0.950 

0.974 

1.009 

0.944 

0.878 

1.000 

All 

24 

1.403 

0.956 

0.951 

0.974 

1.009 

0.944 

0.878 

1.000 


*Indicates significance at a * 0.25 




Table 22: Ratios of Variances About Exponential 

Fit Using Variation of Parameters, Stage 3-4 



Terms Entered 


in Equation 





CO 

rH 

rH 

9k 

o 

00 

1 

so 

«k 

sr 

ON 

H 

rH 

* 

o 

rH 

00 

1 

vO 

<r 

in 

*H 

O 

* 

rH 



rH 

•H 

rH 

H 

rH 

rH 



H 


I unction 

d.f. 

43 

40 

39 

37 

37 

35 

34 

33 

29 


Const. Only 

43 

1.000 









Lin. T , T , D. 
x* m* L 

40 

1.122 

1.000 








Lin. T T D , P 
x m l. r 

39 

1.148 

1.017 

1.000 







Lin. & Int. w/o P 

r 

37 

1.089 

0.971 

0.955 

1.000 






Lin. & Sq. w/o P 
n r 

37 

1.152 

1.026 

1.009 

1.057 

1.000 





Lin. & Sq. w/ P r 

35 

1.130 

1.007 

0.990 

1.037 

0.981 

1.000 




All w/o P f 

34 

1.284* 

1.144 

1.125 

1.179 

1.115 

1.136 

1.000 



Lin. & Int. w/ P r 

33 

1.079 

0.962 

0.946 

0.991 

6.937 

0.955 

0.841 

1.000 


Ail 

29 

1.187 

1.058 

1.040 

1.090 

1.031 

1.050 

0.925 

1.100 

1.000 


indicates significance at a * 0.25 


%% 



Table 23: Ratios of Variances About Exponential 

Fit Using Variation of Parameters, Stage 4-5 


Terms Entered 
in Equation 


m 

m 

m cn eg co m 

r—4 t— 4 «— 4 i— 4 *— 4 »— 4 r— 4 


Function 

d.f. 

30 

30 

29 

29 

29 

28 

27 

27 

26 

25 

24 

24 21 


Const. Only 

30 

1.000 












Lin. Only 

30 

1.000 

1.000 











Lin. T Only 
m 

29 

1.001 

1.001 

i 000 










Lin. T^ Only 

29 

1.208 

1.208 

1.20/ 

1.000 









T , D. , & T D_ 
x* L x L 

29 

1.208 

1.208 

1.207 

1.000 

1.000 








Lin. T^, T | 

28 

1.196 

1.196 

1.194 

0.989 

0.989 

1.000 







Lin. T x . T . D l . P r 

27 

1.153 

1.153 

1.152 

0.954 

0.954 

0.964 

1.000 






Lin. & Int. w/o 

27 

1.184 

J.184 

1.182 

0.979 

0.979 

0.989 

1.027 

1.000 





Lin. & Sq. w/o P r 

26 

1.379* 

1.379 

1.377 

1.141 

1.141 

1.153 

1.196 

1.164 

1.000 




All w/o P 

r 

25 

1.460 

1.460 

1.458 

1.208 

1.208 

1.221 

1.266 

1.232 

1.059 

1.000 



Lin. & Sq. w/ P^. 

24 

1.273 

1.273 

1.271 

1.503 

1.503 

1.065 

1.104 

1.075 

0.923 

0.872 

1.000 


Lin. & Int. w/ P 

r 

24 

1.216 

1.216 

1.214 

1.006 

1.006 

1.017 

1.054 

1.026 

0.882 

0.832 

0.955 

1.000 

All 

21 

1.117 

1.117 

1.116 

0.924 

0.924 

0.934 

0.969 

0.943 

0.810 

0.765 

0.878 

0.919 1.000 


♦Indicates significance at a = 0.25 

Table 2A : Ratios of Variances About Polynomial 

Fit Using Generalized Least Squares, Stage 0-1 
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0k 

r 
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m 

CN 

CN 
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m 


m 

rH 

•k 

0k 

0k 

* 

i 

i 

i 

i 

l 

1 

i 

i 

rH 

rH 

H 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 


Function d.f. 238 237 237 237 235 235 234 232 232 230 229 228 224 


Const. Only 

238 

1.000 












Lin. D l Only 

237 

1.021 

1.000 











Lin. T Only 
m 

237 

1.198* 

1.172 

1.000 










Lin. T x Only 

237 

1.167 

1.142 

0.974 

1.000 









T , D. , & T D. 
x* L* x L 

235 

1.215 

1.190 

1.015 

1.042 

1.000 








Lln - V V °L 

235 

1.225 

1.199 

1.023 

1.050 

1.008 

1.000 







Lln ‘ V V P r 

234 

1.222 

1.196 

1.020 

1.048 

1.005 

0.998 

1.000 






Lin. f» Int. w/o 

232 

1.219 

1.193 

1.018 

1.045 

1.003 

0.995 

0.997 

1.000 





Lin. 6 Sq. w/o P f 

232 

1.396 

1.366 

1.166* 

1.197 

1.148 

1.140 

1.142 

1.145 

1.000 




Lin. & Sq. w/ P f 

230 

1.394 

1.364 

1.163 

1.194 

1.146 

1.1,38 

1.140 

1.143 

0.998 

1.000 



All w/o P r 

229 

1.474 

1.442 

1.231 

1.264 

1.212 

1.203 

1.206 

1.209 

1.056 

1.058 

1.000 


Lin. L Int. w/ P 

r 

228 

1.075 

1.052 

1.879 

0.921 

0.884 

0.877 

0.879 

0.881 

0.770 

0.771 

0.729 

1.000 

All 

224 

1.424 

1.393 

1.189 

1.221 

1.171 

1.163 

1.165 

1.168 

1.020 

1.022 

0.966 

1.325 1.000 


♦Indicates significance at a ■ 0.25 


Table 25: Ratios of Variance About Polynomial 

Fit Using Generalized Least Squares, Stage 1-2 


Terms Entered 
in Equation 
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* c 

•— i i-h m 
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00 

1 
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H 
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H 

•H 

rH 

rH 

rH 

H 

H 

rH rH 


Funct ion 

d.f. 

43 

42 

42 

42 

40 

40 

39 

37 

37 

35 

34 

33 29 


Const. Only 

A3 

1.000 












Lin. Only 

A2 

1.394* 

1.000 











Lin. T Only 
m 

A2 

1.097 

0.751 

1.000 










Lin. T Only 

42 

0.983 

0.706 

0.940 

1.000 









V V 4 T A 

40 

1.302 

0.934 

1.244 

1.324 

1.000 








V V d l 

40 

1.509 

1.083 

1.442 

1.534 

1.159 

1.000 







Lin. T , T , D , P 
x* m* . r 

39 

1.464 

1.050 

1.399 

1.488 

1.124 

0.970 

1.000 






L,in. 6 Int. w/o P^ 

37 

1.950 

1.398* 

1.863 

1.982 

1.498 

1.292 

1.332 

1.000 





Lin. 6 Sq. w/o P 

r 

37 

1.327 

0.952 

1.268 

1.349 

1.019 

0.879 

0.906 

0.681 

1.000 




Lin. 6 Sq. w/ P^ 

35 

1.305 

0.936 

1.247 

1.327 

1.003 

0.865 

0.892 

0.669 

0.984 

1.000 



All w/o P f 

34 

2.269 

1.628 

2.168 

2.307 

1.174 

1.504 

1.550 

1.164 

1.710 

1.738 

1.000 


Lin. f» Int. w / P 

r 

33 

2.697 

1.935 

2.578 

2.743 

2.072 

1.787 

1.842 

1.383* 

2.033 

2.066 

1.189 

1.000 

All 

29 

2,378 

1.705 

2.271 

2.417 

1.826 

1.575 

1.624 

1.219 

1.791 

1.821 

1.047 

0.881 1.000 


♦Indicates significance at a ■ 0.25 


Table 26: Ratios of Variance About Polynomial 

Fit Using Generalized Least Squares, Stage 2-3 
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rH 

rH 

rH 

rH 

rH 

rH 

rH rH rH 


Function 

d.f. 

38 

37 

37 

37 

35 

35 

34 

32 

32 

30 

29 28 24 


Const. Only 

38 

1.000 











Lin. Only 

37 

1.452* 

1.000 










Lin. T Only 
n 

37 

0.980 

0.675 

1.000 









Lin. T^ Only 

37 

0.976 

0.672 

0.995 

1.000 








T , D , 6 T D 
x* L* x L 

35 

1.437 

0.989 

1.466 

1.473 

1.000 







Lin. T , T , D 
x m L 

35 

1.431 

0.985 

1.460 

1.467 

0.996 

1.000 






Lin. T , T„, D. , P 
x* M* L r 

34 

1.434 

0.988 

1.464 

1.471 

0.998 

1.003 

1.000 





Lin. (m Int. v/o 

32 

1.372 

0.945 

1.400 

1.407 

0.955 

0.959 

0.956 

1.000 




I.in. 6 Sq. w/o P f 

32 

1.439 

0.991 

1.468 

1.475 

1 . 001 

1.006 

1.003 

1.048 

1.000 



Lin. 6 Sq. w/ P^ 

30 

1.402 

0.966 

1.431 

1.438 

0.976 

0.980 

0.977 

1.022 

0.975 

1.000 


All w/o P 

r 

29 

1.522 

1.048 

1.553 

1.561 

1.059 

1.064 

1.061 

1.110 

1.058 

1.086 

1.000 

Lin. 6 Int. w/ P^ 

28 

1.324 

0.912 

1.351 

1.358 

0.921 

0.926 

0.923 

0.965 

0.920 

0.944 

0.870 1.000 

All 

24 

0.993 

0.683 

1.013 

1.017 

0.691 

0.693 

0.692 

0.723 

0.690 

0.708 

0.652 0.749 1.000 


*Indicates significance at a • 0.25 


Table 27: Ratios of Variance About Polynomial 

Fit Using Generalized Least Squares, Stage 3-3 


Terms Entered 


In Equation 
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Function 

d.f. 

43 

42 

42 

42 

40 

40 

39 

37 

37 

35 

34 

33 29 


Const. Only 

43 

1.000 












Lin. Only 

42 

0.973 

1.000 











Lin. T Only 
m 

42 

0.996 

1.019 

1.000 










Lin. T Only 

42 

1.150 

1.176 

1.154 

1.000 









V V 4 T A 

40 

1.154 

1.181 

1.159 

1.004 

1.000 








Lin. T x . T m . D l 

40 

1.129 

1.155 

1.133 

0.982 

0.978 

1.000 







Lin. T . T , D P 
x m l r 

39 

1.139 

1.165 

1.143 

0.991 

0.987 

1.009 

1.000 






Lin. 4 Int. w/o P^ 

37 

1.081 

1*. 106 

1.086 

0.941 

0.937 

0.958 

0.949 

1.000 





Lin. ft Sq. w/o P f 

37 

1.100 

1.125 

1.104 

0.956 

0.953 

0.974 

0.966 

1.017 

1.000 




Lin. L Sq. w / P 
1 r 

35 

1.110 

1.135 

1.114 

0.965 

0.962 

0.983 

0.975 

1.026 

1 . 009 

1.000 



Ail w/o P 

r 

34 

1.062 

1.086 

1.066 

0.924 

0.920 

0.941 

0.933 

0.982 

0.966 

0.957 

1.000 


Lin. ft Int. w/ P^ 

33 

1.003 

1.025 

1.006 

0.872 

0.868 

0.888 

0.880 

0.927 

0.911 

0.903 

0.944 

1.000 

All 

29 

1.029 

1.053 

1.033 

0.895 

0.892 

0.912 

0.904 

0.952 

0.936 

0.927 

0.969 

1.027 1.000 


Ul 
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Table 28: Ratios of Variance About Polynomial 

Fit Using Generalized Least Squares, Stage 4-5 
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In Lquation 
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Func 1 1 on 

d. f . 

30~ 

^■4 

30 

7. 2 J 

29 

H 

I2CTI 

H 

2~ 

H 

1 

1 

r-f 

1 

*— « 

i 

«— « 

i 

*-« 

Const. Only 

30 

1 . 000 

J J 

27 

26 

25 

sr 

24 

2l“~ 

Lin. Dj Only 

3u 

1 . 000 

1 . 000 











Lin. T Only 

29 

1.245 

1.245 

1 . 000 










V V ‘ T x D L 

29 

1.245 

1.245 

1 . 000 

1 . 000 









Un - V V D l 

28 

1.230 

1.230 

0.988 

0.988 

1 . 000 








Un - V V b L- f . 

27 

1 .186 

1.186 

0.953 

0.953 

0.965 

1 . 000 







Lin. ft Int. w/o P 

r 

27 

1.211 

1.211 

0.973 

0.973 

0.985 

1.021 

1 . 000 






Lin. It Sq . w/o P 

r 

20 

1. 359* 

1.359 

1.092 

1.092 

1.105 

1.146 

1.121 

1 .000 





All w/o P 

r 

25 

1.615 

1.615 

1.298 

1.298 

1.31 3 

1 . 362 

0.956 

1.189 

1 . 000 




Lin. ft Sq. w / P 

r 

24 

1.158 

1 .158 

0 . 9 30 

0.930 

0.941 

0.976 

1 . 333 

0.852 

0.717 

1 . 000 



I.ln. & Tnt. w/ P 

r 

24 

1.261 

1.261 

1 .013 

1.013 

1 .025 

1.063 

1.041 

0.928 

0. 781 

1.089 

1 .000 


All w/ P 

r 

21 

1.357 

1.357 

1 .090 

1 . 090 

1.103 

1.144 

1 . 120 

0.998 

0.840 

1.172 

1.076 

1.000 


^Indicates significance at • j - 0.25 



Table 29 : Ratios of Variance About Exponential 

f it Usln;» Generalized Least Squares, Stap t e 0—1 


TERMS EnTEREO 
i n Equation 


m 


OJ 


iTs 

l 


cc 

I 

VO 


-=r 

i 


Ov 


irv 

i 


LTn 

I 


i r. 

i 


?2 


229 


228 


2.c4 


C ^ 


9r / 

X- >' 

r — — ' 


Fl'NCTI ON 
Const. Only 

d . r . 

23P 

c y 

1 .000 

, 









u 




Lin. D, Only 

^37 

.972 

1.000 










Lin. T x Only 

237 

1.111* 

1.143 

1 .000 









V °L> V D L 

235 

1 .249 

1 .285 

1 . 124* 

1.000 








Lin. T x * T m , 

235 

1 . 190 

1 .224 

1 .071 

.953 

1 .000 







Lin. T x , T m , Dj_, P r 

234 

1. 192 

1 .226 

1.073 

.995 

1 .002 

1 .000 






Lin. 1 |NT. W/O 

232 

1.030 

1.059 

.927 

.824 

.865 

.864 

1 .000 





Lin. & SQ W/C P R 

232 

i.6o4 

1.650 

1.444 

1.284* 

1 .34? 

1 .346 

1*558 

1.000 




All V# / G P p 

230 

1.605 

1.651 

1 .444 

1.285 

1 .348 

1.346 

1-559 

1 .000 

1 .000 



Lin. £ SQ w/P 

R 

229 

1.478 

1 .520 

1.330 

1 . 1 83 

1 .242 

1 .240 

1.U35 

.921 

.921 

i .000 


Lin. «u Int. w/P r 

228 

1 .232 

1.268 

1 . 109 

.987 

1.035 

1.034 

1.197 

.768 

.768 

.834 

1 .000 

All W/P R 

224 

1-537 

1.58! 

1.383 

1.321 

1 .292 

1 .289 

1.493 

.956 

•95 p 

1 .o4o 

1.247 1.000 


LP 

fj 


♦Indicates significance at a = O.25 


Table 30: Ratios of Variance About Exponential 

Fit using Generalized Least Squares, Stage 1-2 
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1 

1 
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1 

1 
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r— 


I 



Function 

d. f. 

.._s 

42 

42 

40 

40 

3 ? 

37 

37 

35 

_> .. 

?>i 29 

Const. Only 

**3 

1 .000 











Lin. D|_ Only 

42 

1 . 394 * 

1.000 










Lin. T x Only 

42 

.964 

.706 

1.000 









x » °L' T x #i \ 

40 

1.226 

.879 

1.246 

1.000 








LtN * V T m' d l 

40 

1.576 

1.132 

i.6o4 

1.287 

1.000 







Li«. T x , T m , D l P r 

39 

1 .610 

1.155 

1.636 

1.313 

1.020 

1.000 





un 

Lin. &. Int. w /0 Pr 

37 

1.075 

. .771 

1.092 

.876 

.681 

.667 

1.000 





Lin. & SQ W/0 P R 

37 

1.496 

1.073 

1.520 

1 .220 

• 948 

.929 

1*393 

1.000 




All W /0 P r 

35 

1.580 

1.133 

1 .606 

1.289 

1 .001 

.981 

1.470 

1.056 

1 .000 



Lin. & SQ w/P R 

3 ^ 

2.693 

1.931* 

2.737 

2.197 

1.707 

1 .672 

2.505 

1 .800 

1.705 

1.000 


Lin. & Int. W/P r 

33 

14.395 

3.152 

4.466 

3*565 

2.785 

2.729 

4.090 

2.937 

2.782 

1.632* 

1 .000 

All w/P r 

29 

2.676 

. r . - r\ 

1.919 

c*r 

2.720 

2.183 

1 .696 

1 .662 

2.491 

1.769 

1 .694 

,9<jk 

.609 1.000 


Table 31: Ratios of Variance About Exponential 

Fit Using Generalized Least Squares, Stage 2-3 


Terms Entered 
In Equation 
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-=r 


OJ 


o 



rn 


IT'. 
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o 



I 


ITN 


cc 

I 


co 

VO 


Oa 

I 


I 

o 


LT\ 

I 


I 


F UMC T I ON 


o-r. 3ft 


21 


3I 1 35 35 3 1 * 32 32 3^ ?9 2* 2-'« 


Const. Only 

33 

1 .000 




Lin. Only 

37 

1.474* 

1.000 



Lin. T Only 

X 

37 

.976 

.662 

1 .000 


T x > 0^, T x # Dj^ 

35 

1.470 

.997 

1.507 

1.000 

Lin. T , T , D, 
X M L 

35 

1.467 

.995 

1.504 

• 998 

Lin. T , T , D. , P 
X M L R 

3 1 * 

1.473 

1 .000 

1.510 

1.002 

Lin. & Int. w /0 P r 

32 

1 .424 

.966 

i. 46 o 

.969 

Lin. & SQ w /0 P R 

32 

1 . 44 i 

.973 

1.477 

• 960 

All w /0 P 

r 

30 

1.405 

•953 

l. 44 o 

.956 

Lin. & SQ w/P R 

29 

1.566 

1.063 

1.605 

1.065 

Lin. & Int. w/P r 

28 


.911 

1-377 

.914 

All w/P r 

24 

•993 

.674 

1.017 

.675 


•Indicates significance at a =0725 


1.000 







1 .oo 4 

1 .000 






• 971 

.967 

1.000 





.982 

.978 

• 975 

1.000 




• 958 

• 95 1 * 

.987 

•975 

1 .000 



1.067 

1.063 

1.099 

1.087 

1.115 

1.000 


.916 

.912 

.943 

.932 

.956 

.858 

1.000 

.677 

.674 

.697 

.689 

.707 

.634 

.739 1.000 


%% 



Table 32: Ratios or Variance About Exponential 

Fit Using Generalized Least Squares, Stage 3-^ 


Terms Entered 
in Equation 




(VI 

•v 
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10 

1 

CO 

«* 

0 

•x 

1 

00 

1 

VO 

»x 

-=r 
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0 

•s 

00 

1 

vo 

•V 

-=t 

1 

ON 

1 

ITN 

1 

0 

ITN 

1 

ITN 

1 


F UNCT 1 ON 

D.F. 

4 3 

42 

40 

40 

39 

37 

31 

35 

34 

33 

29 


Const. Only 

43 

1.000 











Lin. Dl Only 

42 

.931 

1.000 










Lin. T x Only 

42 

1.145 

I.23O 

1.000 









Tx> T x* D L 

4o 

1.146 

1.231 

1.001 

1 .000 








Lin. T x , T m , D|_ 

4o 

1.121 

1.204 

•979 

.978 

1.000 







L,n# T x* V D L> P r 

39 

1.137 

1.222 

•99 1 * 

.992 

1.015 

1.000 






Lin. &. Int. w /0 P 

37 

i.o4i 

1.118 

.910 

.908 

.929 

• 915 

1.000 





Lin. &. SQ w/O P R 

37 

1.142 

1.227 

.998 

.997 

1.019 

1.005 

1.097 

1.000 




All w/O P r 

35 

1.118 

1.200 

.976 

•975 

•997 

•983 

1.074 

•978 

1.000 



Lin. & SQ w/P 

R 

34 

1.211 

1.301 

1.058 

1.056 

1.080 

1.065 

1.163 

1.060 

1.083 

1.000 


Lin. & Int. w/P r 

33 

.860 

.92** 

.751 

.750 

.767 

.756 

.826 

•753 

.769 

,v° 

1.000 

All w/P 

' R 

29 

1.008 

1.O83 

.881 

.880 

.900 

.887 

.969 

.883 

.902 

.832 

1.172 1.000 


Table 33 : Ratios of Variance About Exponential 

Fit Using Generalized Least Squares, Stage 4-5. 


56 


each phenological interval regardless of its environment. It is against this 
assumption that the subsequent functional forms were tested. Inspection of 
the tables supports several general conclusions: 

1. ) fits in the plant-emerge interval were significantly bettered 

only when the nonlinear effects were included. Note, how- 
ever, that a linear equation using temperature only was 
"near significant." Assuming that the variance will not 
change significantly, inclusion of more data points to in- 
crease degrees of freedom could easily make this case 
significant; 

2. ) fits in the emerge-joint interval were significant when 

linear temperature functions were used. Further significant 
improvement was not noted until squared terms were added to 
the function; 

3. ) The joint-head interval was highly sensitive to linear day- 

length functiouo . The addition of interaction terms signi- 
ficantly bettered Lhe fit. This was the only interval In 
which the addition of precipitation further bettered the fit 
significantly; 

4. ) the linear daylength terms were significant in the head- 

soft dough interval. No significant differences were noted 
in this interval by adding terms to this linear function; 

5. ) interestingly, for the most part the soft dough-ripe interval 

showed no significance due to the addition of any environ- 


mental variables; 
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6.) in almost no instances were the fits significantly bettered 
by the inclusion of precipitation in the models. 

Based upon these facts, and the conclusions reached earlier concerning 
the relative merits of the exponential and polynomial models, it would 
appear that the following functional forms merit further, more detailed 
study. 


Stage 0-1: R ■ H + H~T 

1 2 x 


R - H. + H~T + H_T + H,D T + H,T 2 + H,T “ + H Q D T 2 
1 2x 3m 4 L 6x 7m 8L 


Eq. 34 


Stage 1-2: R - H, + H_T 

1 2 x 


2 2 2 

R * H + H 0 T + H-T + H,D t + H,T + H,T + H R D 
1 2x 3m 4 L 6x ' m ° L 


Eq. 35 


Stage 2-3: R =» H, + H.D T 

1 4 L 


R - H. + H.T + H~T + H,D. + H. n T D T + H..T D, + H. _T T 
1 2x 3m 4L 10 x L 11 m L 13 x m 

R * H. + H,T + H_T + H. D + HJ* + H in T D_ 

I 2x 3m 4L 5r 10 x L 

+ H..T D. + H 10 P D. +H n TT + H..T P„ + H...T P„ 

II m L 12 r L 13 x m 14 x r 15 m r 


Eq. 36 


Stage 3-4: R ■ H. + H.D T 

1 4 L 


Eq. 37 


Stage 4-5: R ■ 


Eq. 38 


R ■ H- + H 0 T + H.D, + H. n T D T 
1 2 x 4L 10 x L 


Tests on Independent Data 

The tests described in the preceeding section, using variance ratios and 
F-tests, can be used only to draw conclusions concerning the least squares 
fit for the data used in these fits. However, the data used in the least 
squares procedures represents, by necessity, averages over time and 
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position. In order to assess the appropriateness of the resulting para- 
meters in a predictive mode, these parameters must be applied to independent 
data not included in the least squares procedures. 

This was done on a limited scale. In order to test in this manner, 
daily envimomental data must be prepared and averaged spatially. Since this 
is a very time consuming process, time limitations for this contract phase 
precluded the inclusion of many location/years of data in this process. 

Three location/years of data were prepared for each phenological interval for 
th'',e tests. A simple testing program was written to apply the parameters 
obtained from the least squares procedures to these daily environmental 
v riables, and cumulative stage computed independently within each interval. 

A listing of this test program, together with typical output, is included in 
Appendix D. 

Even though the sparcity of data prepared for this test program allowed 
no meaningful statistical tests, some useful insights were gained. In all 
instances, when higher order terms were included in the equations, the 
resulting predicted stage was highly erratic. The reasons for this are 
fairly obvious. The parameters are estimated using temporally averaged data. 
Thus, the least squares is applied to data for which extreme values have 
been "damped out." When these parameters are applied to daily data, many of 
the values will be outside the range used to estimate the parameters, and 
nonlinear effects will appear exaggerated. When the linear functions were 
tested, these erratic results did not appear. Figures III-IX show plots for 
the location/years tested in each stage interval using linear functions. Of 
particular interest are the plots for the second interval, representing 



Figure Ills Tests of Linear Function, Stage 0-1 
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emergence to jointing. Overall, the success of the modeling in predicting was 
quite good. Of further interest is the reaction of the model to the original 
dormancy assumptions. It appears that the model can indeed predict a 
dormancy period, indicated by falling growth rates, and still adequately 
predict the jointing date. 

Based upon the tests described in this section, it is recommended that 
a linear function of the form 

R = H. + H 0 1 + H^T + H,D. Eq. 39 

1 2 x 3 m 4 L 

be utilized for further testing in the LACIE program. The parameters 
suggested for further testing are shown in Table 34. 


Table 34: Polynomial Coefficients Suggested 


Stage 


H i 


H 2 

H 3 

H 4 



0-1 

-0.014919 

0.0038970 

0 . 

0 . 

1-2 

-0.0003y918 

0.00043509 

0 . 

0 . 

2-3 

-0.216419 

0 . 

0 . 

0.019021 

3-4 

0.314583 

0 . 

0 . 

-0.018610 

4-5 

0.244711 

0.0046211 

0.0015439 

-0.022684 


SUMMARY , CONCLUSIONS AND RECOMMENDATIONS 
Summary 

Winter wheat phenological information W3S gathered for crop reporting 
districts located in the Great Plains and Rocky Mountain regions of the 
United States. Corresponding meteorological data were obtained from the 
National Climatological Center for these location/years, for values of daily 
maximum temperature, minimum temperature, and precipitation. Daylength 


value! were computed. 
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These data, after editing and reduction, were utilized in least squares 
fitting programs using both the method of variation of parameters, and the 
method of simultaneous adjustment of observations and parameters (generalized 
least squares). Parameters were estimated using these techniques for 
various forms of polynomial, exponential, and Robertson's "triquadrat ic" 
functions. The functions resulting from these fits were then tested for 
goodness of fit using a- posteriori variances from these fits. Further, the 
resulting functions were tested in a predictive mode by applying them to 
limited data sets withheld from the least squares procedures. 

C onclusions 

The investigations conducted support of the following general 
conclusions : 

1. ) the use of phenological reports at the crop reporting 

district appears viable. The use of such information re- 
quires no separate data gathering, as it is accomplished 
within the U.S.D.A. reporting service. The data, however, 
are not without disadvantages. Few states report all 
phenological stages, and many report other than the 
"standard" stages (plant, emerge, joint, head, soft dough, 
ripe). No states, to this investigator's knowledge, report 
information concerning dormancy or spring greenup, 

2. ) because only end points are availanle, at best, in the 

phenological reports, only average growth rates may be 
computed and utilized for parameter estimation within each 
•stage interval for each location/year. This, in turn, 
necessitates that the corresponding environmental values 


be averaged spatially and temporally over the CR1) and the 
stage interval; 

4. ) tor the emerge to joint interval special considerations are 

necessary. This interval Includes the long winter dormancy 
period, and some criteria must be established in data pre- 
paration which reflects this fact. Later testing procedures 
indicated that the behavior of the resulting models were 
relatively insensitive to minor variations in the dormancy 
assumpt ions made ; 

5. ) the generalized least squares modeling techniques, combined 

with numerical methods for linearization and t met ion 
evaluation, appear to provide a powerful and flexible tool 
for parameter estimation within a variety of function types. 
While some difficulties were encountered in obtaining con- 
vergence for models with higher order terms, due to the 
approximat ions involved, the problems were relatively minor, 
and are certainly reconcilable. Fairly good parameter ap- 
proximations are necessary in using the model. Further the 
method lacks efficiency in testing of variance for inclusion 
of terms when compared with stepwise regression. 

However, the flexibility of the technique, including 
its ability to accommodate nonlinear models, and various 
configurations within each model is a great advantage. 

Further the fact that the method recognizes and includes, 
computationally the statistical variability within all 
observed quantities, and minimizes with respect to all these. 


Is folt to result In more realistic parameter estimates than 
those resulting from regression analysis; 
b.) the "triquadratic" model, advocated hv Robertson (O tor spring 
wheat, shows promise for winter wheat. The model is para- 
metrically efficient, and may easily accommodate nonlinear and 
interaction effects. However, the model is highly nonlinear 
in the parameters, and requires special handling in order to 
separate the davlength and temperature terras, necessary in 
order to avoid singular normal equation matrices. Further, 
the model requires good parameter approxim.it ions , or risks 
converging to unrealistic parameter sets. These problems were 
not completely solved during this phase oi the contract. 

7.) for the more straight forward functional forms, it was found 
that an exponential form ottered few advantages over the 
simpler and more easily applied polynomial forms; 

8.1 In virtually no instance was the inclusion of precipitat ion 
within l he mode l f ound t o be s i gn i t lean t ; 

9.) tests of functions on Independent data revealed that the 

more highly nonlinear functional forms produced very erratic 
results when driven with dailv environmental variables. This 
was not unexpected, since parameter estimates were generated 
using averaged data, and the dally values used for testing 
often lay outside the range of these averaged values, resulting 
in greatly magnified results for those points; 


r 
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Recommendations 

The following recommendations appear to be viable for the immediate 
future based upon the results of this study. 

1. ) as of this writing, the following functions are recommended 

for inclusion and testing within the LACIE program: 

Stage 0-1: R - -0.014919 + 0.0038970 T x 

Stage 1-2: R * -0.00039918 + 0.00043509 T x 

Stage 2-3: R * -0.216419 + 0.019021 D 

Eq. 40 

Stage 3-4: R * 0.314583 - 0.018610 D L 

Stage 4-5: R * 0.244711 + 0.0046211 T + 0.0015439 T 

x m 

-0.022684 D 

2. ) as an obvious first step, additional location/years of data 

should be prepared and included in the least squares modeling 
program. These data should be carefully chosen to represent 
geographic positions and crop years not already included. 

These efforts are underway; 

3. ) many more location/years of data should be prepared and 

added to the testing program, in order that sufficient re- 
dundancy will result to permit statistical testing from the 
results of this program. This necessitates an extensive 
data-gathering and editing effort, but is of extreme 
importance, and is underway; 

4. ) investigations should be conducted concerning ways to drive 

the model other than with daily environmental data values. 

It is suspected that if cumulative or moving averages of 
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environmental variables are used, the more highly nonlinear 
models will not react so erratically, and may offer a viable 
alternative to the linear functions recommended in 1.) above; 

5. ) Further testing should be conducted concerning the cumulative 

effects of errors within the model. A test program should 
be written which drives the model throughout the entire crop 
year for each location, and outputs error within each stage 
as well as cumulative error through each stage; 

6. ) Work should proceed on the development of a "triquadratic" 

model suitable for winter wheat; 

7. ) variance propagation techniques should be utilized in order 

to assist in the daily use of the crop calendar by analyst 
interpreters within the LACIE program. The end results of 
such studies would, ideally, yield a date at which the crop 
reaches each phenological stage and an error interval about 
that date. It is important in pursuing this goal to recognize 
that the variances of the parameters, as well as those of the 
independent variables, be included in such studies. Since 
the output from the generalized least squares procedures 
results in a full variance/covariance matrix for the para- 
meters, this can be accomplished. 
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